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The dynamical cluster-decay model (DCM) is developed further for the decay of hot and rotating 
compound nuclei (CN) formed in light heavy-ion reactions. The model is worked out in terms of 
only one parameter, namely the neck-length parameter, which is related to the total kinetic energy 
TKE(T) or effective Q- value Qe//(T) at temperature T of the hot CN, defined in terms of the 
CN binding energy and ground-state binding energies of the emitted fragments. The emission of 
both the light-particles (LP), with A<4, Z<2, as well as the complex intermediate mass fragments 
(IMF), with 4<A<20, Z>2, is considered as the dynamical collective mass motion of preformed 
clusters through the barrier. Within the same dynamical model treatment, the LPs are shown 
to have different characteristics as compared to the IMFs. The systematic variation of the LP 
emission cross section aLP, and IMF emission cross section oimf, calculated on the present DCM 
match exactly the statistical fission model predictions. It is for the first time that a non-statistical 
dynamical description is developed for the emission of light-particles from the hot and rotating CN. 
The model is applied to the decay of ^^Ni* formed in the ''^S-|-^'*Mg reaction at two incident energies 
£c.m. =51. 6 and 60.5 MeV. Both the IMFs and average TKE spectra are found to compare resonably 
nicely with the experimental data, favoring asymmetric mass distributions. The LPs emission cross 
section is shown to depend strongly on the type of emitted particles and their multiplicities. 

PACS numbers: 25.70.Jj, 23.70.-|-j, 24.10.-i, 23.60.-|-e 



I. INTRODUCTION 

Light compound nuclei (CN) formed in low energy 
(E/A<15 MeV/nucleon) heavy-ion reactions are highly 
excited and carry large angular momenta. The com- 
pound systems so formed decay by emitting multiple light 
particles (n, p, a) and 7-rays, which in the statistical 
Hauser-Feshbach (HF) analysis is understood as an emis- 
sion process from the equilibrated CN, resulting in CN 
fusion cross sections. Apparently, the decay process must 
depend on temperature and angular momentum depen- 
dent potential barriers 0. For light compound systems 
with AcN > 40, the above noted light-particles (LP) 
emission is always accompanied by intermediate mass 
fragments, the IMFs of Z>2 and 4<A<20, also called 
"complex fragments" or "clusters", whose contribution, 
though small of the order of 5 to 10%, is to be included in 
the CN fusion cross sections. Then, the temperature- and 
spin-dependent potentials must also be mass-asymmetry 
dependent. In other words, the structure effects of the 
compound system also become important. 

In order to understand the IMF production, not only 
the HF analysis is extended to include the fragments 
heavier than a-particle in the BUSCO code Q or in 
the Extended Hauser-Feshbach scission-point model [3, 
but also other statistical fission model descriptions ^ 
have been used that are based on either the scission- 
point or saddle-point configuration, in the GEMINI code 

tor the saddle-point "transition-state" model (TSM) 
B H S 0, respectively. The LP emission in these 



fission models is still treated within the statistical HF 
method. It is interesting to mention that light-nucleus 
emission can also be qualitatively well described in the 
framework of a generalized rotating liquid-drop model, 
proposed recently by Royer and collaborators fiol ITU . 
Since the measured angular distributions and energy 
spectra of emitted complex fragments are consistent with 
fission-like decays of the respective compound system, 
the fusion-fission process is now well established in light 
dinuclcar systems 0. The statistical fission models, 
for which the fission-decay of a CN is determined by 
the phase space (level density) available at the saddle- 
point 0,0 (or scission-point (3|) configuration, may how- 
ever lack in terms of not including more explicitly the 
structure effects of the compound system. Large struc- 
ture effects have been shown to be important in the 
^^Ni coi npound system through a strong resonance be- 
haviour 12, 13] of the excitation functions of large-ande 
^*Si+^*Si clastic and inelastic scattering yields |S, 14|. 
Although neither similar resonant effects nor orbiting 
processes 1151 have been evidenced in the ^^S-|-^''Mg re- 
action [idllTllTsIl studied here, a fully dynamical theory 
for more complete description of the emission of both the 
LPs and IMFs in the framework of the statistical model 
of the decay of such a hot and rotating nuclear ^^Ni sys- 
tem remains highly desirable. 

A recent attempt towards the dynamical treatment of 
the decay of a hot and rota ting nucleus is the work of 
Gupta and collaborators [11 Hfl [H IH mHll , where a 
dynamical collective clusterization process is proposed as 
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a possible alternative of the fission process. Both the LPs 
and IMFs are considered as the dynamical mass motion 
of preformed fragments or clusters through the barrier. 
Note that, in terms of the barrier picture, a cluster-decay 
process is in fact a fission process with structure effects 
of the CN also included via the preformation of the frag- 
ments, but without any phase space arguments (i.e. with 
no level density calculations). Alternatively, the dynam- 
ical fission process has been considered by some authors 
123, J26| si mply as a continuous deformation of the CN 
I27L 123 . 129I I30I l3l| . The above noted dynamical cluster- 
decay model (DCM) of Gupta et al. jll S El El is 
so-far used to calculate the decay constant and total ki- 
netic energies (TKEs) of IMFs alone, and that too in a 
re-normalization procedure only. In other words, only 
the angular momentum £=0 (s-wave) solution and TKE 
for fixed £-value are considered. In this paper, we develop 
the DCM further for the calculation of actual (summed 
up, "total" £) cross sections and the average total ki- 
netic energies TKEs for both the emitted light particles 
and complex IMFs. Interesting enough, the light parti- 
cles, though treated within the same dynamical collective 
clusterisation process, are found to possess different char- 
acteristic properties. Brief reports of the pres ent work 
have recently been presented elsewhere |2^l24j. 

The data for IMFs chosen here, to apply the DCM, 
are those of the ^^S-l-^'^Mg^^^Ni* reaction, where the 
mass spectra for A = 12 to 28 fragments, and the aver- 
age total kinetic energy (TKE) for only the most favored 
(enhanced yields) a-nucleus fragments, are measured at 
two incident energies Eiab = 121.1 and 141.8 MeV, or 
equivalently at Ec.m. = 51.6 and 60.5 MeV, respectively 
1^, _16J (only even-even fragments are observed for Ec.m. 
= 51.6 MeV). In this experiment, by detecting both frag- 
ments in coincidence, it was possible to deduce the pri- 
mary mass distribution for the decay process i.e. the 
mass distribution before the occurrence of the secondary 
light-particles emission from the fragments. These pri- 
mary, pre-evaporation mass distributions show that the 
mass-asymmetric channels are favored over the symmet- 
ric ones, with a-nucleus, A=4n, fragments having en- 
hanced yields. The IMF emission cross section, estimated 
in one of these experiments is 59 ± 12 mb. The CN 
fusion cross section data due to multiple LP emission at 
these energies are also deduced later from the same ex- 
periment and as quoted in j7|, are 1080 ± 130 and 
1050 ± 100 mb, respectively, at Ec.m. — 51.6 and 60.5 
MeV, which fit the other available earlier measurements 
at similar and other energies for this system [HIHIll. 
The total fusion cross section is then the sum of this cross 
section due to the LP emission and the fission-like IMF 
emission cross section. 

The paper is organised as follows. Section II gives 
briefly a description of the DCM for hot and rotating 
light-mass nuclei. The model is presented for the emis- 
sion of both the LPs and IMFs. The calculations are pre- 
sented in section III. The discussion of our results and a 
summuary constitute section IV. 



II. THE DYNAMICAL CLUSTER-DECAY 
MODEL FOR HOT AND ROTATING 
COMPOUND SYSTEMS 

The DCM for hot and rotating nuclei is a reformulation 
of the preformed cluster model (PCM) of Gupta and col- 
laborators for ground-state decays in cluster radioactiv- 
ity (CR) and related phenomena 34, 35, 36, 37, 38, 39j. 
Thus, like PCM, DCM is also based on the dynamical (or 
quantum mechanical) fragmentation theory of cold fusion 
phenomenon in heavy ion reactions and fission dynani- 
^ M M yi including the prediction of CR 
43I lil ] . This theory is worked out in terms of: 
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(i) the collective coordinate of mass (and charge) asym- 
metry (H and L stand, respectively, for heavy and light 
fragments) 

i]={Ah-Al)/{Ah + Al) 



(and Tjz = (Zh - Zl)/{Zh + Zl)), 

and (ii) relative separation R, 

which in DCM characterizes, respectively, 

(i) the nucleon-division (or -exchange) between outgoing 
fragments, and 

(ii) the transfer of kinetic energy of incident channel 
{Ec.m.) to internal excitation (total excitation or total 
kinetic energy, TXE or TKE) of the outgoing channel, 
since the fixed R = Ra (defined later), at which the pro- 
cess is calculated, depends on temperature T as well as 
on 77, i.e. i?(T, ry). This energy transfer process follows 
the relation (see Fig. 1) 



Ec 



QoutiT) 



TKE{T) + TXE{T). 

(1) 

Note that for the 32s+24Mg^56]s^j* taction, Q„ is posi- 
tive (=16.68 MeV) and hence adds to the entrance chan- 
nel kinetic energy Ec.m. of the two incoming nuclei in 
their ground states, and Qout{T) is negative and differ- 
ent for different exit channels at a fixed temperature T. 
The CN excitation energy E^j^ and its temperature T 
(in MeV) are related as 



E*c^^{A/%)T^-T. 



(2) 



Using the decoupled approximation to R- and ry- 
motions, the DCM defines the decay cross section, in 
terms of partial waves, as [2lll2^ 



1=0 



2fiEc. 
If 



(3) 

where, Pq, the preformation probability, refers to rj- 
motion and P, the penetrability, to R-motion. Appar- 
ently, for £=0 (s-wave) ao = ^PqP, which is an equiv- 
alent of decay constant A = i^oPoP (or decay half-life 
Ti/2 = In2/X) with i/q as the barrier assault frequency. 
In other words, tJo and A differ through a constant only. 
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Thus, like in PCM, here the complex fragments (both 
LPs and IMFs) are treated as the dynamical collective 
mass motion of preformed clusters or fragments through 
the barrier. The structure information of the CN en- 
ters the model via the preformation probabilities Pq (also 
known as the spectroscopic factors) of the fragments. 

The Pq is given by the solution of stationary 
Schrodinger equation in 77, at a fixed R, 



{- 



d 1 d 



dr] dT] 



+ V{R,r],T)}r{v) = E''r{v), 

(4) 

with ^=0,1,2,3... referring to ground-state (i/ = 0) and 
excited-states solutions. Then, the probability 



Po(AO -I ^(7?(A0) 1^ 



'A' 



(5) 



where i= 1 or 2, the heavy (H) and light (L) fragment, 
respectively, and for a Boltzmann-likc function 



(6) 



v=0 



The constant R = Ra = Ct = Ci + C2 in fixed 
empirically as the first turning point of the penetration 
path (like in Fig. 1) for the ground state {g.s, T=0) 
decay, since this value of R (instead of CN radius Rq) 
assimilates the effects of the deformations of two frag- 
ments and neck formation between them [S] . Thus, for 
the deformation effects included, we have for g.s. decay 

2 

Ra{T^Q)^Ct or, in general, ^Ct + ^dR{P,), 

1=1 

(7) 

an equivalent of lowering of the barrier for deformed 
fragments This is important here because the ex- 

pected shapes of some of the observed fragments in the 
exit channel of the reactions studied here could be de- 
formed [141 li^ . Thus, the deformation effects of the 
fragments (and the neck formation between them), in- 
cluded here within the extended model of Gupta and 
collaborators are treated via a neck-length 

parameter ^ SR{(3i) at the scission configuration which 
simulates the two centre nuclear shape parameteriza- 
tion, used for both the light ^3 Ell and heavy 0| 
nuclear systems. This method of introducing a neck- 
length parameter y]SR is similar to that used in both 
the scission-point "Sl and saddle-point 0, statistical 
fission models. The alternate, i.e., the calculation of frag- 
mentation potential V^rj) and scattering potential V(R) 
for deformed nuclei, thou gh shown to be difficult [45«| . 
is being worked out [i^ Isoj. Also, we use here the 
Siismann central radii Ci ~ Ri ~ b^/Ri (in fm), where 
Ri = 1.28Ay^ - 0.76 + 0.8A;^^^ fm and surface thick- 
ness parameter b =0.99 fm. Note that the Ct are different 
for different 7y- values and hence Ra{T = 0) also depends 
on rj. 



For the decay of a hot CN, we use the postulate |2C 
I21I |2^ for the first turning point 



RaiT) = Raiv, T = 0) + Ai?(77, T), 



(8) 



depending on the total kinetic energy TKE(T). The cor- 
responding potential V{Ra) acts like an effective Q-value, 
Qeff, for the decay of the hot CN at temperature T, to 
two exit-channel fragments observed in g.s. T=0, defined 
by 

QeffiT) = B(T)-[Bl{T = Q) + Bh{T = 0)] 

= TKE{T)^V{Ra{T)). (9) 

Here B's are the respective binding energies. Thus, AR 
in Eq.® gives the change in TKE(T) w.r.t. TKE(T=0), 
which is taken to depend on 77. As a first approximation 
to Eq. @ , in our earlier calculations 21, 23,] , we used 
a constant average AR, which also takes care of the addi- 
tional J2 SRiPi) effects of the deformations of fragments 
and neck formation between them, 



RaiT) = Ctiv,T) + ARiT). 



(10) 



Note that here Ct is also taken to depend on temperature, 
to be defined later. In the following calculations, we shall 
first use a constant AR{T) and then show the effect of 
using the actual 7/-dependent AR(ri,T), calculated from 
Qeff or V{Ra,T) defined by Eq.®. 

The above defined decay of a hot CN to two cold (T=0) 
fragments, via Eq.® , could apparently be achieved only 
by emitting some light particle(s), like n, p, a, or 7-rays 
of energy 



E,. = B{T)-BiO) = QeffiT)-Qout{T = 0) 
= TKE{T) - TKE{T = 0) 



(11) 



which is zero for the g.s. decay, like for exotic CR. 
Note that the second equality in Eq. is not de- 

fined for a negative Qout{T — 0) system since the nega- 
tive TKE(T=0) has no meaning. Apparently, Eq. l(TT)) 
w.r.t ® suggests that the emission of light-particles 
starts early in the decay process. The exit channel frag- 
ments in ® are then obtained in the ground-state with 



TKE(T=0), as can be seen by calculating E\ 



CN 



E^, 



E*c^^E., =1 Qout{T) I +TKE{T = 0)+TXE{T). (12) 

The excitation energy TXE(T) is used in, not treated 
here, the secondary emission of light particles from the 
fragments which are otherwise in their ground states with 
TKE(T=0) in the radial motion. Instead, we compare 
present calculations with the primary pre-secondary- 
evaporation fragment emission data |y, lla |. Thus, by 
defining Qeff{T) as in ®, in this model we treat the LP 
emission at par with the IMF emission. In other words, 
a non-statistical dynamical treatment is attempted here 
for the first time for not only the emission of IMFs but 
also of multiple LPs, understood so-far only as the statis- 
tically evaporated particles in a CN emission. It may be 
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reminded here that the statistical model (CN emission) 
interpretation of IMFs is not as good as it is for the LP 
productionH H H IE II El 

The collective fragmentation potential V{R,ri,T) in 
Eq.Q) is calculated according to the Strutinsky method 
by using the T-dependent liquid drop model energy 
Vldm of Ref. IsJ, with its constants at T=0 re-fitted 
|2ll |23 | to give the recent experimental binding ener gies 
[S^ . and the "empirical" shell corrections SU of Ref. |53 
(In the Appendix of '2l'| and Eq. (8) of ^2^, 0^=0.5, in- 
stead of unity). Then, including the T-depcndence also 
in Coulomb, nuclear proximity, and ^-dependent poten- 
tial in complete sticking limit of moment of inertia, we 
get 



V{R,r,,T) = Y,[Vldm{A,Z,,T) 

exp{ 



i=l 



i=l 



Tn 



+ E,{T) + Vp{T) + Vt{T). (13) 

The shell corrections SU{T) in p3|l are taken to go to 
zero exponentially with T, with To = 1.5 MeV |li|. The 
other T-dependent terms Ec, Vp and Vi are defined as 
follows: 



E,{T) = ZHZLeyR{T), 



(14) 



where the charges Zi are fixed by minimizing the po- 
tential V{R,ri,T) in the charge asymmetry coordinate 
T]z (for fixed R and i] coordinates). Note, here the T- 
dependence of Ec is the same as in Ref. |51|- Next, 



Vp{T) = A7TR{T)jb{T)^s{T)), 



(15) 



with R{T) and $(s(r)), respectively, as the inverse of 
the root mean square radius of the Gaussian curvature 
and the universal function, independent of the geometry 
of the system, given by |53 



R{T) 



Ci{T)C2{T) 
Ct{T) ' 



(16) 



^, _ , "^(s - 2.54)2 - 0.0852(s - 2.54)^ s < 1.2511 
)) - i ^3 437ga.p(__s_) ^ > 1.2511 



(17) 



and 7, the nuclear surface energy constant, given by 

2' 



7 = 0.9517 



1 - 1.7826 



N- Z 
A 



MeVfm-^. (18) 



The overlap or separation distance, in units of b, between 
the two colliding surfaces 



siT) 



R-Ct{T) 
b{T) 



with the temperature dependence of radii Ri , used in Ci , 
taken from Ref. |^ as 



i?,(r) = 1.07(1 + 0.01T)A| 
and surface width [56| . 

h{T) = 0.99(H-0.009T2). 



(19) 



(20) 



Here the proximity potential Vp of Ref. |53 is extended 
to include the T-dependencc. Similarly, 



(21) 



2/(r) ' 

where, in the complete sticking limit, 

/(T) = /s(T) = fiRl + IahtuCI + \AL'mGl. (22) 

Here [i = [AhA^I^Ah + AL)]m = \Am{\ - 77^) is the 
reduced mass, with m as the nucleon mass. Note that 
this limit is defined for the separation distance Ai?, or 
Ai?, to be within the range of nuclear proximity (<2 fm). 

Finally, the fm(ja;-value in Eq. Q is the critical £-value, 
in terms of the bombarding energy Ec.m., the reduced 
mass ^ and the first turning point Ra of the entrance 
channel r]in , given by 



4 - RaVMEc.m. -v{Ra,mn,e ^ o)]/h, 



(23) 



or, alternatively, it could be fixed for the vanishing of 
fusion barrier of the incoming channel, called if us or 
else the £-value where the light-particle cross section 
aLP 0. This, however, could also be taken as a vari- 
able parameter 0, 0| . 

The mass parameters 5^^(77), representing the kinetic 
energy part in Eg. Q) , are the smooth classical hydrody- 
namical masses [57|, since we are dealing here with tem- 
peratures where the shell effects are almost completely 
washed out. 

The P in Eq. © is the WKB integral. 



P = exp[~- j ^ {2ti[ViR) - Q,ff]Y/^dR], 



(24) 



solved analytically ,35|, with Rh as the second turning 
point (see Fig. 1) satisfying 



V{Ra) = V{Rb) ^Qeff = TKE{T). 



(25) 



Note that, since we do not know how to add the ^-effects 
in binding energies, the ^-dependence of Ra is defined by 



ViRa) = QeffiTJ=0), 



(26) 



i.e. Ra is the same for all ^-values, given by the above 
equation, and that Qeff{T,£) = V{Ra,i)- Then, using 
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(|25|l . Rb{i) is given by the ^-dependent scattering poten- 
tials, at fixed T, 

V{R, T, e) - ZHZLe^/R{T) + Vp{T) + 1/,(T), (27) 

which is normalised to the exit channel binding en- 
ergy. This means that all energies are measured w.r.t. 
Bl{T) + Bh{T). Such a potential is illustrated in Fig. 1 
for the decay of ^^m* at T=3.39 MeV {E^.m.=^l-& 
MeV), using different ^-values. The first turning point 
Ra is marked for R — Ra — Ct + Ai?-values and the 
second turning point Rb for only the £ = case of 
Ra — Ct + Ai?. Note that as the £-value increases, the 
(5e//(=TKE)-value increases, since the decay path for all 
the i?- values begins ai R = Ra- For the reaction under 
study, at Ec.m. =51.6 MeV, the critical angular momen- 
tum £c = 32?i from Eq. (|23|l and, for the vanishing of 
barrier for incoming channel, if us = 49/i. It is interest- 
ing to note that this £c value is only slightly lower than 
the value of icrit = 34fi. that has been extracted I p, | a 
from the complete fusion cross section data [T^ l32L l33l| 
using the sharp cutoff approximation. The correspond- 
ing values for the higher energy i?c.m. =60. 5 MeV (T=3.60 
MeV) case will be given later, at the appropriate place. 

III. CALCULATIONS 

The reaction 32S-|-24Mg^56]v^j* ^las been studied at 
two incident energies [Ec.m. =51.6 and 60.5 MeV) which 
correspond to 1.8 and 2.2 times the Coulomb barrier 
^, 16] . As already mentioned in the Introduction, at 
such energies, the fission-like IMF cross section <jimf 
is about 6% of the LP emission cross section a^p. 
Specifically, oimfIolp = (59 ± 12)m6/(1050 ± 50)m6 
at Ec.m. =60.8 MeV In the following, however, 

we choose to use the data of Ref. IQj, and the above 
noted energy of Ref. [3 is close to one of the energies 
of the chosen data. In the experiment of Ref. 0, for 
the fission-like decay of ^^Ni*, a complete mass spec- 
trum of IMFs is measured, beginning with mass Al — 12 
fragment. Later on the measurements of the IMF mass 
spectrum have been extended at Eiab = 130 MeV (equiv- 
alently, Ec.m. = 55.7 MeV) in order to include the ^Be 
fragment yields [s^- In this experiment [s^, an en- 
hanced emission yield for ^Be by a factor of 1.5 to 1.8, 
over the two a-particles, is observed, which in the Ex- 
tended Hauser-Feshbach method (EHFM) is shown 
to be related to an increased deformation of the heav- 
ier fragment *^Cr. More recently, ternary events from a 
conjectured hyperdeformed ^^Ni CN have been observed 
at Eiab = 165 MeV {equivalentlpEcm. = 70.7 MeV) 
[5^. The two sets of data (Ref. 0^ and ^5^) are con- 
sistent with the strong deformation effects found in ^^Ni, 
as is also popula ted by the ^^Si+^*Si fusion-evaporation 
reaction |60ll61|. 

Figure 2 gives for ^^Ni* the mass fragmentation po- 
tentials V{A) at different ^-values, for the fixed T=3.39 
MeV and i? = i?a = Ct 1.28 fm. The R- value is chosen 



for the best fit to the IMF cross section data (see below) . 
Two interesting results are apparent in this graph: (i) 
Because of T-dependent Vldm, the non-a, Z=N even-A, 
fragments also appear at minima which are in addition to 
the preferred a-nucleus structure present in the macro- 
scopic liquid drop energy due to the " Wigner term" . The 
shell corrections 5U are already nearly zero at these tem- 
peratures. In other words, with the addition of tempera- 
ture in the potential, not only the shell structure effects 
vanish but also there is no explicitly preferred a-nucleus 
structure left. It is worth noting that the same behav- 
ior is also known for fission calculations based upon 
either the saddle-point picture 0, 0, or the scission- 
point picture 0- (ii) The structure in the fragmentation 
potential (the positions of minima and maxima) is in- 
dependent of the £-value, though very important effects 
of symmetric or asymmetric mass division are present 
here in Fig. 2. Apparently, the favored (lower in en- 
ergy) asymmetric mass distribution at zero and smaller 
^-values go over to the symmetric one for partial waves 
with angular momenta i near the ^c-value. In particular, 
at lower ^-values the light particles (plus the correspond- 
ing heavy fragments) are strongly favored, over the heav- 
ier fragments (IMFs), but this situation gets reversed at 
or near the ^c-value. As we shall see in the following, 
these results have important consequences for the rela- 
tive contributions of LP and IMF emission at different 
ts to the total decay cross section. It may be relevant 
to mention here that we already know from the experi- 
ments of Beckerman et al. |H3 that the emission of IMFs 
starts only beyond a certain energy (and hence beyond a 
certain angular momentum) and that, for lower energies, 
only LP emission occurs which gives the complete fusion 
cross section. This result is already given by the DCM 
worked out in s-wave (^=0) approximation |20l I21L |23| . 

Figure 3 shows the calculated preformation factors 
Po{£,Al), the penetration probability P(i,AL) and the 
cross section ae^A^), with £ summed over imax = ^c, for 
use of the fragmentation potentials of Fig. 2. Only the 
energetically most favored mass fragments are considered 
for both the LPs and IMFs. First of all, two important 
results can be drawn for Po{Al) (shown as solid thin 
line in Fig. 3): (i) Po{Al) is a strongly oscillating func- 
tion with maxima only at Z=N, even-A fragments; (ii) 
the preformation yields are large for light-particles and 
the asymmetric fragments. In other words, in agreement 
with experiments JJJ , an asymmetric mass distribu- 
tion of IMFs is favored by the preformation factors. On 
the other hand, the penetrability P{Al) (dashed line) is 
almost a monotonically decreasing function, with P being 
relatively small for symmetric fragments. Thus, P{Al) 
also support the asymmetric mass distribution. The to- 
tal cross section (t{Al) = Y^e^o ^ii-^L) for each fragment 
(solid thick line ) is a combined effect of these two terms, 
to be discussed below in comparison with experiments. 
Here we notice that its behavior is given more by Pq{Al) 
than by P{Al). In the following, we first study the vari- 
ation of these quantities with angular momentum £. 
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Figure 4 shows the variation of Pq with £ for the ener- 
getically favored LPs {A < 4) and the even-A, N=Z IMFs 
(the contribution of the energetically unfavored odd-A 
IMFs is small at all fs). The maximum ^- value, £max, 
is taken to be eqaul to the ^c-value. We notice that, 
whereas Pq decreases for LPs with an increase of £, it 
increases for IMFs as £ increases and then starts to de- 
crease at a large value. Po{£) for "^He behaves like the 
LPs and that the behaviour of all LPs (A<4) is different 
from that of the IMFs (A>4). Also, for heavier IMFs 
(A>16), Pq is almost zero for £ < 18h. Furthermore, in 
Fig. 5, the P's for LPs are also large, rather the largest 
{P — 1 for proton emission since there is no barrier at all 
^'s), but the same for £ < 18fi is nearly zero for all IMFs. 
Thus, for the penetrability P also, the behavior of LPs is 
different from those of the IMFs. This result for the cross 
sections means that the lower £ values contribute mostly 
to the LP cross section aLP and that the higher £ values 
{£ > 18h) contribute to the fission-like IMFs production 
cross section ajMF- This is illustrated in Fig. 6, where 
CTiP, (JiMF and the total cross section 

C^Total = CTLP + CflMF 

are plotted for each £. (The values of the three cross 
sections for £max = £c are also given in the brackets of 
the legend). We notice that ctlp is already zero for £ > 
31h (a value close to £c), which means that not only the 
two processes of LP and IMF emissions get separated 
at £ « 18h, but also the decay process stops aX £ — tc- 
value (at least for the LPs in the present calculations). 
The IMF cross section starts to contribute only for t > 
18S and is maximum bX £ — £(.. Furthermore, if the £ 
coordinate is extended up to, say, £fus (not shown in the 
figure), the contribution to LP emission remains zero but 
the IMF cross section goes on increasing as £ increases. 
This means that if the decay process continues beyond 
£ ~ £c, only <timf contributes to axotai- In other words, 
in DCM, £max = £c seems to be an automatic choice, 
fixed by the initial conditions of the experiment, as in 
Eq. H23|) . Alternatively one could choose £max at an 
i?- value where cr^p — s- 0, as is discussed below. 

The individual contributions of IMFs are illustrated in 
Fig. 7, for the ^Be emission and the "^^Si emission. We 
notice that the contributions of the lighter IMFs towards 
(J IMF are much larger than that for heavier IMFs. This 
result is consistent with the observation of an asymmet- 
ric mass distribution, which is favored. It is interesting 
to note that the same results (as presented in Figs. 6 and 
7) can be obtained in the statistical fission model calcu- 
lations for this reaction (see, for example. Fig. 14 in |^). 
The noticeable difference is that, in the statistical fission 
model 0,^3) (^fission (= o'imf) also reduces to zero at 
£ = £c, which is due to the chosen phase space (the sharp 
cutoff approximation) in that model. Another point of 
interest to note in DCM is that the so-called promptly 
emitted LPs are really not that prompt but they do have 
a considerable overlap with the binary-decay process (of 
cluster emission) for the higher £-values. This is also con- 



sistent with the statistical fission model of |a, [ig • It is 

true that LP emission starts early but continues along 
with the emission of IMFs till the decay process itself 
stops for £ > £c. 

Before giving the quantitative comparison of ctlp and 
<^iAi F calculations with experimental data, we study here 
the role of changing Ai?- value and the nature of emitted 
light-particle(s). Fig. 8 shows the results of DCM with 
the use of different Ai?-values for the LPs and IMFs, 
taking Ai? = 0.41 fm for LPs, but keeping the same 
A^=1.28 fm (as in Fig. 6) for IMFs. We notice that the 
magnitude of (total ^-summed) a^p reduces considerably 
(by a factor of about 2) where as ct/m f remains nearly the 
same (rather an increase by about 20% can be observed) . 
Furthermore, if we also change the proton-emission to 
neutron-emission, as in Fig. 9, the ^c-value remains the 
same, but the magnitude of ctlp reduces further to about 
60%, keeping the gimf almost unchanged. These results 
in Figs. 8 and 9 are to be compared with the respec- 
tive measured values of (Tlp=1050 mb and (T/j\/p=60 mb. 
The agreement for uimf can be further improved if the 
^-values are summed only upto the point where ulp 0; 
then the calculated cr/MF = 103 or 106 mb in Fig. 8 or 
Fig. 9, respectively. Note that the drastic reduction of 
(Tlp in Fig. 9 occurs because the lower ^-values {t < 3h) 
also do not contribute to the LP emission cross section. 
It is known for experiments that it is more difficult to 
evaporate a neutron than a proton, as is also shown to 
be energetically the case in Fig. 2. However, we do not 
make a search for the exact emission of LPs, since we do 
not attempt a one-to-one comparison with the data. On 
the other hand, the general success of DCM is demon- 
strated here for the first time to give the LP emission in 
a non-statistical formalism. 

The total cross section a = J^'^i for decay of ^^Ni* 
at Ec.m. =51.6 MeV to light fragments Al (for both the 
LPs and even-A, N=Z IMFs) is plotted in Fig. 10 for 
the case of Fig. 8, and compared with other calculations 
based either on the saddle-point transition-state model 
(TSM) 6] or on the scission-point model, the so-called 
EHFM scission-point model 3] , and available experimen- 
tal data mill for > 12 IMFs. The TSM calculations 
for LPs are performed within the HF formalism 6], and 
hence are the same for EHFM model IS]. The EHFM 
calculations for IMFs are done for Al >12 fragments 
and are thus joined straight from Al =4 to Al =12. 
Only even-A fragments are plotted since the IMF-spectra 
at Ec.m. =51.6 MeV is measured for only even-A frag- 
ments 0, ^3 . For the LP emission at this energy, as al- 
ready stated above, the measured (fusion or evaporation 
residue) cross section is available (ctlp = 1050 ± 100 mb) 
but the separate yields for each emitted LP is not given 
for a possible direct comparison between the experiment 
and model calculations. 

In Fig. 10, we can notice that the known disconti- 
nuity at the point between Al =4 and 6 in both TSM 
and EHFM calculations, due to the use of HF formal- 
ism for LPs {A <4), is no more present in the present 
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DCM calculations. The DCM treats both the LPs and 
IMFs emission in a similar manner, although the present 
calculations (corresponding to the case of Fig. 8) over- 
estimate glp by a factor of more than two. This dis- 
crepancy in DCM is mainly due to the contribution of 
mass-one-particle, here the proton. For example, the re- 
placement of proton with neutron (Fig. 9) yields the 
calculated glp to be in better agreement with the data. 
Apparently, for the calculation of the evaporation residue 
cross section u^p or a^p + crjMF, it is very important to 
know exactly the contributing particles (i.e. their multi- 
plicities) for comparisons with experiments. Note further 
that the HF analysis gives nearly equal cross sections for 
each of the four emitted particles {Al — 1-4:), whereas a 
decreasing function of the light-particle mass is obtained 
in DCM. It will be of high interest to measure the LP's 
trends in the near future. 

For the IMFs, in Fig. 10, the general comparison be- 
tween the experimental data and DCM is of the same 
quality as for the TSM or EHFM, at least for Al < 22. 
For Al > 22 the TSM and EHFM predictions appear 
to be better. Note that in DCM calculations A^=0.41 
fm for LPs and Ai?=1.28 fm for IMFs. However, a 
closer comparison of DCM calculations with experiments 
(Fig. 12, discussed below) favors the use of a fragment- 
dependent AR or the actual AR{ri), calculated from 
Qeff, presented in Fig. 11 for light mass fragments. 

In Fig. 11, we notice that AR{ri) has an oscillatory 
nature, if compared to the smooth variation of Ct or 
Ct -\- AR with rj. The maxima in AR{ri) correspond 
to a-nuclei IMFs and the minima to N=Z non-a-nuclei 
IMFs, with the odd-A^ fragments lying in between. For 
light-particles, the AR{r]) values increase almost mono- 
tonically. In any case, the division between the LPs and 
the IMFs is clearly evident from the variation of AR with 
rj. 

Figure 12 shows the DCM calculations for use of dif- 
ferent average Ai?- values and the actual AR{ri) obtained 
in Fig. 11 from calculated Qeff- These calculations are 
presented here only for even- A IMFs. We have also added 
here the DCM calculations for AS=1.28 (from Fig. 10). 
It is clear that Ai?=1.28 fm gives the optimum fit to 
IMF data, though the oscillatory nature of data is al- 
most smoothed out, particularly for the heavier IMFs. 
Interesting enough, this oscillatory structure of the cross 
section gets restored with the use of actual AR{ri) ob- 
tained from calculated Qeff, though the fit with data is 
now deteriorated. Apparently, an improvement in AR{ri) 
is required. This calls for an improvement in the calcu- 
lations of Qeff and hence in the ground-state and T- 
dependent binding energies. 

The histograms in Fig. 13 show the comparisons of the 
absolute IMF cross sections for the best fit (Ai?=1.28), 
the experimental data 0, and the two alternate 
model-calculations of TSM (|| and EHFM 0], for Al > 
12 at both the available energies. Similarly, Fig. 14 
shows the DCM calculated excitation functions (cross 
sections at different Ec.m.) for the emission of ^^C from 



the excited ^^Ni* CN. We notice in Fig. 14 that, in- 
dependent of the choice of Ai?-value (i.e. a constant or 
T-dependendent value), the ^^C emission cross section, 
cr(^^C), increases as the incident energy increases and 
reaches a maximum around Ec.m. — 90 MeV and then 
starts to decrease at higher incident energies. It is in- 
teresting to note that similar results are obtained in the 
HF calculations, using BUSCO code, for the emission of 
i^C from "4-118^^* jgj. From Figs. 13 and 14 it is 
clear that the DCM contains the required features of the 
experimental data, as well as of other models (EHFM Q 
and TSM 6]). A better treatment of the binding en- 
ergies and missing aspects, such as the deformations of 
the fragments and neck formation between them, would 
also ameliorate the predictions of the present model. In 
view of this hope, in the following, we further analyze the 
comparisons of the DCM calculations for average TKEa 
with the experimental data of Ref. 6] . 

Figure 15 shows the DCM calculated average total ki- 
netic energy 

TKEiAL) - ^^j^TKEii,AL), 

compared with the experimental data for the ^^S-|-^^Mg 
reaction leading to hot ^^Ni* at the two chosen ener- 
gies. Here, for each fragment, the TKE for each £ is av- 
eraged over its corresponding production cross section ai 
w.r.t. the total cross section cf{Al) = Y^i^q o-e{AL)- We 
have also calculated the total kinetic energy TKE{Al) 
for a best fit to the £-value. Apparently, the comparisons 
with data are reasonably good for both the calculations 
and it is difficult to distinguish between the calculated 
TKE and TKE. The maximum i'-value is nearly the 
same in both case. However, it is not clear why this max- 
imum £- value is much less than the value. The simple 
model-dependence used for handling the deformations of 
the fragments and neck formation between them need 
further improvements. 

IV. SUMMARY 

We have further developed the DCM for the decay of 
a hot and rotating CN, formed in light heavy-ion reac- 
tions, into multiple LP evaporation and IMF emission. 
The LP emission (evaporation residue) cross section aLP, 
constitutes the CN fusion cross section (J fusion for a neg- 
ligible emission of IMFs, since afusion = o^lp + (^imf, 
also referred to as UTotai- The statistical equilibrated 
CN evaporation process, successfuU for the emission of 
LPs, could not explain the IMF emission. Alternatively, 
the IMF emission alone could be understood as the sta- 
tistical fissionprocess in the saddle-point |^ or scission- 
point model (si. On the other hand, in DCM, both the 
LPs and IMFs are treated identically as the dynamical 
collective mass motion of preformed fragments or clusters 
through a barrier, i.e. quantum mechanical tunnelling of 
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clusters that are considered pre-born with different prob- 
abihties before they actuaUy penetrate the barrier. Thus, 
the cluster preformation probabilities contain the struc- 
ture effects of the CN, that are found to be important in 
the description of the measured excitation functions of 
large-angle elastic and inelastic scattering yields in the 
experiments under study. 

The DCM is worked out in terms of only one parame- 
ter, the neck-length parameter, that depends on the to- 
tal kinetic energy of the fragments TKE(T) at the given 
temperature T of the CN, which itself is defined for the 
first time in terms of the binding energies of the emitted 
fragments in their ground-states and the binding energy 
of the hot CN. The hot CN is considered to achieve its 
ground-state by giving away its extra binding energy to 
the emitted LPs, which is shown to leave the emitted 
IMFs in their respective ground states with total kinetic 
energy TKE(T=0). The remaining (excitation) energy 
TXE(T) must go in to the emission of secondary light- 
particles from the IMFs which are otherwise already in 
their ground-states in the radial motion. Such an emis- 
sion of secondary light-particles is not included here in 
the DCM so-far; rather the model predictions are com- 
pared with the primary IMF experimental data |^, cor- 
rected for such an emission. 

The DCM is apphed here to the decay of ^^Ni*, formed 
in the '^^S-t-^'^Mg reaction at two incident energies Ec.m.— 
51.6 and 60.5 MeV, where both the LP cross section 
and IMF spectra, as well as the total average kinetic 
energy (TKE) for only the favored a-nucleus fragments, 
are measured j^. The interesting result of DCM is that 
both the preformation factors and penetrabilities, as a 
function of angular momentum, behave differently for the 
LPs and the IMFs. In other words, there is an explicit di- 
vision at mass-four fragment between the LPs and IMFs 
with *He belonging clearly to the LP regime. The prefor- 
mation factor is shown to contribute more to the observed 



behaviour of IMF cross section crjMF, which can be com- 
pared with the experimental data reasonably well, favor- 
ing an asymmetric distribution. Furthermore, the vari- 
ation of both a LP and (Jimf with angular momentum, 
as well as the individual contributions of IMFs to aiMF, 
and the excitation functions of the emitted IMFs, match 
exactly the predictions of the statistical fission model, 
and the HF analysis. Since, unlike fission models, the 
DCM does not depend on the chosen phase space it has 
the advantage that the ^mas-value is fixed by the initial 
conditions of the experiment via £c, rather than by the 
available phase space. This distinguishing feature is ev- 
ident in a IMF not going to zero when cr^p goes to zero 
at tmax = (-c- The comparison of (Tlp, however, depends 
strongly on the type of particles involved and their mul- 
tiplicities, as expected. The calculated TKE also repro- 
duces the experimental data, though at an ^ < £c-value, 
which has perhaps to do with the way the deformations of 
the fragments are included here simply through the same 
neck-length parameter that accounts for the temperature 
effects. The model is being improved both for the ne- 
glected deformation effects and neck-formation between 
them as well as the binding energies used to calculate 
this neck-length parameter. 
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Figure Captions: 

Fig. 1 The temperature and angular momentum 
dependent scattering potentials, illustrated for 
56jvi* _>12 c Ti at T = 3.39 MeV (equiva- 
lently, E^.m. — 51.6 MeV). The potential for each t 
is calculated by using y(i?,T,^) = Ec{T)+Vp{T) + 
V^(r), normalized to exit channel T-dependent 
binding energies Bl{T) + Bh{T), each defined as 
B{T) = Vldm{T) + 5U{T). The decay path, de- 
fined by V{RaJ-) — Qeff{T,£) for each £, is shown 
to begin at Ra — Ct + AR for £ = case, where 
AR is the average over 77 of the different neck- 
length parameters AR{= Ra ~ Ct) calculated for 
V{Ra) = Qeff{T,£ = 0) for all possible fragmenta- 
tions. The critical angular momentum £c- value is 
determined from Eq. H23I) . For other details, see 
text. 

Fig. 2 The fragmentation potentials V(A) for the 
compound system ^^Ni* formed in the reaction 
32S-|-24Mg at Ec.m. =51.6 MeV, calculated for dif- 
ferent ^-values at a fixed T=3.39 MeV (correspond- 
ing to Ec.m. =51.6 MeV) and Ra ^ Ct + 1.28 fm. 
The fragmentation potential V{A) ~ Vldm{T) + 
5U{T) + E,{T) + Vp{T) + Vt{T) for fixed R = Ra. 
The fc-value is as per Eq. (f^ . 

Fig. 3 The fragment preformation factor Po(^,^l), the 
penetrability P{£,Al) and the decay cross section 
ae{AL), with £ summed over imax — £c = 32h cal- 
culated on Eq. (123 , for the decay of ^^^i* to var- 
ious complex fragments (both LPs and IMFs), us- 
ing the fragmentation potentials in Fig. 2 based on 
DCM. 

Fig. 4 The variation of Pq with £, for both the 
LPs (dashed hues) and even-A, N=Z IMFs (soHd 
lines), calculated on DCM for the compound sys- 
tem ^^Ni*, using the fragmentation potentials in 
Fig. 2. For proton, the calculated Pq values are 
three times of what are shown in the figure. 

Fig. 5 The variation of P with £, for both the 
LPs (dashed hues) and even-A, N=Z IMFs (soHd 
lines), calculated on DCM for the compound sys- 
tem ^^Ni*, using the scattering potentials as in Fig. 
1. For proton, there is no barrier at any £- value and, 
hence P=l for proton. 

Fig. 6 The variation of evaporation residue cross sec- 
tions due to the LPs cr^p (dotted line), the IMFs 
c^iMF (dashed line), and their sum dTotai (solid 
line) with £, for values up to ^c, calculated on 
DCM for the compound system ^^Ni* formed in 
32s+24j^jg reaction at E^.m. = 51.6 MeV. The pa- 
rameter AR (=1.28 fm) is the same for both the 
LPs and IMFs. The cross sections given in the 
brackets are obtained by summing for £max — £c — 
32h. Here the LPs consist of ^'"^H and ^'-^He. 

Fig. 7 The variation of Po(^,^l), P{£,Al) and ^^(Al) 
with £, for the emission of ^Be and ^^5*1 fragments 



from ^^Ni* formed in ^^S-h^^'Mg reaction at E^.m. = 
51.6 MeV, calculated on DCM, as in Fig. 6. For 
^^S'i-decay, the calculated P is ten times and, Pq 
one tenth, of the plotted values. 

Fig. 8 The same as for Fig. 6, but for use of different 
AP-values for LPs and IMFs. 

Fig. 9 The same as for Fig. 8, but for the proton- 
emission replaced by neutron-emission. 

Fig. 10 The same as for Fig. 3, but for only the decay 
cross section (j{Al) and calculated for the case of 
Fig. 8. The DCM calculations are compared with 
two other model calculations EHFM and TSM of 
Refs. ^1 an d 0, respectively, and the experimental 
data for Al > 12 IMFs. 

Fig. 11 The variation of the first turning point Ra with 
light fragment mass A^ , for cases of Ra = Ct, Ra ^ 
C't + AR{=1.28 fm), and the actually calculated Ra 
from V{Ra) = Qeff{T), £=0, for the decay of^^Ni 
at T=3.39 MeV. Note that the actually calculated 
Ra from V{Ra) — Qeff{T), 1=0 can be written as 
Ra = Ct + AR{ri) where AR{ri) is found to be, in 
general, positive. For some light fragments {Al=1- 
5,8), we use £ >0 since calculated Qeff{T) is larger 
than the barrier for £=0 case. 

Fig. 12 The same as for Fig. 10, but^ DCM 
alone, calculated for different average Ai?-values 
and the actual AR{Al) determined from V{Ra) = 
Qeff{T,£ = 0) (Fig. 11). The DCM calculations 
are compared with the experimental data taken 
fromRef. 0,13. 

Fig. 13 The histograms of the calculated IMF cross 
sections a(Ai)on DCM, compared with the exper- 
imental data and two other model calcula- 
tions EHFM 3 and TSM 0, for even-A^ IMFs 
at Ec.m. =51.6 MeV and for both the odd- and 
even-Ai IMFs at Ec.m. =60.5 MeV. The predicted 
Al =14 fragment cross section on DCM is very 
large, ten times of what is plotted here. 

Fig. 14 The DCM excitation functions i.e. the cross sec- 
tions at different incident cm. energies for emission 
of ^^C fragment from the excited ^^Ni* compound 
nucleus, calculated for a constant and an arbitrary 
T-dependent Ai?-value. 

Fig. 15 The measured 6] and DCM calculated aver- 
age total kinetic energy (TKE) for the reaction 
^^S-l-^^Mg^^^Ni* Al+Ah, at two incident en- 
ergies Ec.m. =51.6 and 60.5 MeV. Also, the total 
kinetic energy TKE for the best fit to ^-value is 
plotted. The average Ai?=1.28 and 1.29 fm, re- 
spectively, for the two energies. 
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